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Abstract. Let A be a minor-closed class of labelled graphs, and let Q n be a random graph 
sampled uniformly from the set of n-vertex graphs of A. When n is large, what is the 
probability that Q n is connected? How many components does it have? How large is its 
biggest component? Thanks to the work of McDiarmid and his collaborators, these questions 
are now solved when all excluded minors are 2-connected. 

Using exact enumeration, we study a collection of classes A excluding non- 2-connected mi- 
nors, and show that their asymptotic behaviour may be rather different from the 2-connected 
case. This behaviour largely depends on the nature of dominant singularity of the generating 
function C[z) that counts connected graphs of A. We classify our examples accordingly, thus 
taking a first step towards a classification of minor-closed classes of graphs. Furthermore, we 
investigate a parameter that has not received any attention in this context yet: the size of the 
root component. It follows non-gaussian limit laws (beta and gamma), and clearly deserves a 
systematic investigation. 



1. Introduction 

We consider simple graphs on the vertex set {1, . . . , n}. A set of graphs is a class if it is 
closed under isomorphisms. A class of graphs A is minor-closed if any minoiQof a graph of A is 
in A. To each such class one can associate its set £ of excluded minors: an (unlabelled) graph 
is excluded if its labelled versions do not belong to A, but the labelled versions of each of its 
proper minors belong to A. A remarkable result of Robertson and Seymour states that £ is 
always finite [32] ■ We say that the graphs of A avoid the graphs of £. We refer to [6] for a study 
of the possible growth rates of minor-closed classes. 

For a minor-closed class A, we study the asymptotic properties of a random graph Q n taken 
uniformly in A n , the set of graphs of A having n vertices: what is the probability p n that Q n is 
connected? More generally, what is the number N n of connected components? What is the size 
S n of the root component, that is, the component containing 1? Or the size L n of the largest 
component? 

Thanks to the work of McDiarmid and his collaborators, a lot is known if all excluded graphs 
are 2-connected: then p n converges to a positive constant (at least 1/^/e), N n converges in law 
to a Poisson distribution, n — S n and n — L n converge in law to a discrete distribution. Details 
are given in Section [3j 

If some excluded minors are not 2-connected, the properties of Q n may be rather different 
(imagine we exclude the one edge graph...). This paper takes a preliminary step towards a 
classification of the possible behaviours by presenting an organized catalogue of examples. 

For each class A that we study, we first determine the generating functions C(z) and A{z) that 
count connected and general graphs of A, respectively. The minors that we exclude are always 
connectecj^] which implies that A is decomposable in the sense of Kolchin |25j : a graph belongs to 
A if and only if all its connected components belong to A. This implies that A(z) = exp(C(z)). 
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^obtained by contracting or deleting some edges, removing some isolated vertices and discarding loops and 
multiple edges 

2 We refer to 28 for an example where this is not the case. 
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Figure 1. A zoo of graphs. Top: the 3-star, the triangle Ks, the bowtie and 
the diamond. Bottom: A caterpillar and the 4-spoon. 

We then derive asymptotic results from the values of these series. They are illustrated throughout 
the paper by pictures of large random graphs, generated using Boltzmann samplers |18| . Under a 
Boltzmann distribution, two graphs of A having the same size always have the same probability. 
The most difficult class we study is that of class avoiding the bowtie (shown in Figure [IJ . 

Our results make extensive use of the techniques of Flajolet and Sedgewick's book |19) : sym- 
bolic combinatorics, singularity analysis, saddle point method, and their application to the 
derivation of limit laws. We recall a few basic principles in Section [2] We also need and prove 
two general results of independent interest related to the saddle point method or, more precisely, 
to Hayman admissibility (Theorems [f7] and [l8) . 

Our results are summarized in Table [T] A first principle seems to emerge: 

the more rapidly C(z) diverges at its radius of convergence p, the more compo- 
nents Q n has, and the smaller they are. 

In particular, when C(p) converges, then the properties of Q n are qualitatively the same as in 
the 2-connected case (for which C(p) always converges [27]), except that the limit of p n can be 
arbitrarily small. When C(p) diverges, a whole variety of behaviours can be observed, depending 
on the nature of the singularity of C(z) at p: the probability p n tends always to 0, but at various 
speeds; the number N n of components goes to infinity at various speeds (but is invariably 
gaussian after normalization) ; the size S n of the root component and the size L n of the largest 
component follow, after normalization, non-gaussian limit laws: for instance, a Gamma or Beta 
law for S n , and for L n a Gumbel law or the first component of a Poisson-Dirichlet distribution. 
Cases where C(z) converges, or diverges at most logarithmically, are addressed using singularity 
analysis (Sections [4] and [5]), while those in which C(z) diverges faster (in practise, with an 
algebraic singularity) are addressed with the saddle point method (Sections [7] to |l0| ) . Section |6] 
gathers general results on the saddle point method and Hayman admissibility. 

Let us conclude with a few words on the size of the root component. It appears that this 
parameter, which can be defined for all exponential families of objects, has not been studied 
systematically yet, and follows interesting (i.e.. non-gaussian!) continuous limit laws, after nor- 
malization. In an independent paper we perform such a systematic study, in the spirit of 
what Bell et al. [4 or Gourdon [22 did for the number of components and the largest compo- 
nent, respectively. This project is also reminiscent of the study of the 2-connected component 
containing the root vertex in a planar map, which also leads to non-gaussian continuous limit 
laws [3]. 

2. Generatingfunctionology for graphs 

Let E be a finite set of (unlabelled) connected graphs that forms an antichain for the minor 
order (this means that no graph of £ is a minor of another one). Let A be the set of labelled 
graphs that do not contain any element of £ as a minor. We denote by A n the subset of A 
formed of graphs having n vertices (or size n) and by a n the cardinality of A n . The associated 
exponential generating function is A(z) = J2 n >o a « z " A 1 " We use similar notation (c„ and C(z)) 
for the subset C of A consisting of (non-empty) connected graphs. Since the excluded minors 
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axe connected, A is decomposable, and 

A(z) =exp(C(z)). 

Several refinements of this series are of interest, for instance the generating function that keeps 
track of the number of (connected) components as well: 

where \G\ is the size of G and c(G) the number of its components. Of course, 

A(z, u) — exp(uC(z)). 
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Table 1. Summary of the results: for each quantity N n , S n and L n , we give 
an equivalent of the expected value (up to a multiplicative constant, except in 
the last line where constants are exact) and a description (name or density) of 
the limit law. The examples are ordered according to the speed of divergence 
of C{z) near its radius p. Spoons are defined in Figure [l] As we get lower in 
the table, the graphs have more components, of a smaller size. The symbol 
PD ( - 1 ' ) (l/4) stands for the first component of a Poisson-Dirichlet distribution of 
parameter 1/4. The results of Appendix B are of slightly different nature since 
they deal with coloured forests. Therefore they do not appear in this table, even 
though they are similar to those of the first two lines. 
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We denote by Q n a uniform random graph of A n , and by N n the number of its components. 
Clearly, 

m _ _ ^M*)' m 

i\[ z n] A (zy w 

and the factorial moment of N n is 

_ [z n }C(zYA(z) 



E(iV„(iV n - 1) • • • (N n - i + 1)) = 



[z n ]A(z) [z n ]A(z 



Several general results provide a limit law for N n if C(z) satisfies certain conditions: for instance 
the results of Bell et al. (¥| that require C(z) to converge at its radius; or the exp-log schema 
of [Ti?l Prop. IX. 14, p. 670], which requires C(z) to diverge with a logarithmic singularity (see 
also the closely related results of [2] on logarithmic structures). We use these results when 
applicable, and prove a new result of this type, based on Drmota et a/.'s notion of extended 
Hayman admissibility, which applies when C(z) diverges with an algebraic singularity. We 
believe it to be of independent interest (Theorem 18 1. 

We also study the size c\ of the root component, which is the component containing the 
vertex 1. We define accordingly 



A{z,v)= ^2 >' 



JG\-i 

ci(G)-l z 



(|G| - 1)!' 

The choice of \G\ — 1 instead of |G| simplifies slightly some calculations. Note that A(z, 1) = 
A'(z) — C (z)A(z). Denoting by S n the size of the root component in Q n , we have 

tw/g l\ c fe a «-fc(fc-D kc k a„_ fc n\ 

V(S n = k) = = -TT7 mT — ■ ( 2 

a n n k\ [n — k) \ a n 

Equivalently, the series A(z, v) is given by 

A{z,v) = C'(zv)A(z). (3) 
The factorial moment of S n — 1 is 

Surprisingly, this parameter has not been been studied before. Our examples give rise to non- 
gaussian limit laws (Beta or Gamma, cf. Propositions 14 or 25 1. In fact, the form ^ of the 



generating function shows that this parameter is bound to give rise to interesting limit laws, 
as both the location and nature of the singularity change as v moves to 1 — e to 1 + e. Using 
the terminology of Flajolet and Sedgewick |19[ Sec. IX. 11], a phase transition occurs. We are 
currently working on a systematic study of this parameter in exponential structures |11| . 

Finally, we denote by C^(z) the generating function of connected graphs of A of size less 
than k: 

CM(*) = 5>^, 

n— 1 

and study, for some classes of graphs, the size L n of the largest component. We have 



(5) 
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3. Classes defined by 2-connected excluded minors 

We assume in this section that at least one minor is excluded, and that all excluded minors 
are 2-connected. This includes the classes of forests, series-parallel graphs, outer-planar graphs, 
planar graphs... Many results are known in this case. We recall briefly some of them, and state 
a new (but easy) result dealing with the size of the root component. The general picture is that 
the class A shares many properties with the class of forests. 

Proposition 1 (The number of graphs — when excluded minors are 2-connected). The 

generating functions C(z) and A(z) = e c ^ are finite at their (positive) radius of convergence p. 
Moreover, the sequence (a n /n\) n is smooth, meaning that na n -\/a n tends to p as n grows. 

The probability that Q n is connected tends to \/A{p), which is clearly in (0,1). In fact, this 
limit is also larger than or equal to l/\/e. The latter value is reached when A is the class of 
forests. 

The fact that p is positive is due to Norine et al. [31J, and holds for any proper minor-closed 
class. The next results are due to McDiarmid [27] (see also the earlier papers [29 , 30 ) . The fact 
that 1/A(p) > l/\/e, or equivalently, that C(p) < 1/2, was conjectured in [30], and then proved 
independently in [1 and [23] . 




FIGURE 2. A random forest of size n = 1165 having two connected components. 



Example 2. A basic, but important example is that of forests, illustrated in Figure [2] We have 
in this case 

where T(z) — ze T ^ counts rooted trees (see for instance |19l p. 132]). The series T. C and 
A = e c have radius p — 1/e, with the following singular expansions at this point: 

T{z) = l -V2(l-ze) 1 /2 +|(1 - ze ) -^(1-ze) 3 / 2 +0((1 - ze) 2 ), 

C{z) = \ ~(l-*e) +2^1(1 -z e )3/2 +0((l- ze ) 2 ), (6) 

A{z) = ^fe -^(1-ze) +^fe^ (1 - zef' 2 +0((1 - zef). 
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The singularity analysis of |19l Chap. VI] applies: the three series are analytic in a A-domain, 
and their coefficients satisfy 

and a n ~ \fec(n). 



Note. When all excluded minors are 2-connected, C{p) always converges, but the nature of the 
singularity of C(z) at p depends on the class: it is for instance (1 — z/p) 3 / 2 for forests (and more 
generally, for subcritical minor-closed classes [15]), but (1 — z/p) 5 / 2 for planar graphs. We refer 
to [H| for a more detailed discussion that applies to classes that exclude 3-connected minors. 

Proposition 3 (Number of components — when excluded minors are 2-connected). 

The mean of N n satisfies: 

E(N n ) ~ 1 + C{fi) 

and the random variable N n — 1 converges in law to a Poisson distribution of parameter C(p). 
That is, as n — > oo, 

W'^ + ^SSr ( 7 ) 

We refer to |271 Cor. 1.6] for a proof. The largest component is known to contain almost all 
vertices, and it is not hard to prove that the same holds for the root component. In fact, the 
tails of the random variables S n and L n are related by the following simple result. 

Lemma 4. For any class of graphs A, and k < n/2, 

P(S„ =n-k) = ^—^ F(L n = n-k). 

Proof. Let us denote by B n the (lexicographically first) biggest component of Q n . Its size is thus 
L„. We have, for n > 2k, 

P(S n =n-k) = ¥(S n = n - k and 1 € B n ) + P(S„ = n - k and 1 £ B n ) 

= ¥{L„ = n - k and 1 € B n ) + P(5„ =n-k and 1 £ B n ) 

= P(l e B n \L n =n- k)¥(L n = n - k) + P(S n = n - k and 1 g B n ) 

n — k „ . 

= F(L n = n-k). 

n 

Indeed, there cannot be two components of size n — k or more. This implies that ¥(S n — 
n — k and 1 ^ B n ) = 0. ■ 

Proposition 5 (The root component and the largest component — when excluded 
minors are 2-connected). The random variables n—S n and n—L n both converge to a discrete 
limit distribution X given by 

'<* = *> = ^ir- 

Proof. By Lemma [4] the two statements are equivalent. The L n result has been proved by 
McDiarmid [57] Cor. 1.6]. 

We give an independent proof (of the S n result), as we will recycle its ingredients later for 
certain classes of graphs defined by non- 2-connected excluded minors. Let k > be fixed. By |2|, 

m , a ,« c„_ fc a fe (" fc 1 ) a k c„_ fc (n-l)!a„_fe 

P(5 n =n- k) = = — — -— - . 

a„ fc! a„_ fc (n-fe-l)!a„ 

By Proposition [l] the term Cn^k/o-n-k, which is the probability that a graph of size n — k is 
connected, converges to 1/A{p). Moreover, the sequence a n /n\ is smooth, so that fer^—fji^ 
converges to /A The result follows. a 
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A more precise result is actually available. Let us call fragment the union of the components 
that differ from the biggest component B n . Then McDiarmid describes the limit law of the 
fragment, not only of his size [271 Thm. 1.5]. 

4. When trees dominate: C(z) converges at p 

Let A be a decomposable class of graphs (for instance, a class defined by a excluding connected 
minors), satisfying the following conditions: 

(1) A contains all trees, 

(2) the generating function D{z) that counts the connected graphs of A that are not trees 
has radius of convergence (strictly) larger than 1 /e (which is the radius of trees) . 

We then say that A is dominated by trees. Some examples are presented below. In this case, the 
properties that hold for forests (Section |3| still hold, except that the probability c n /a n that Q n 
is connected tends to a limit that is now at most 1/y/e. We will see that this limit can become 
arbitrarily small. 

Proposition 6 (The number of graphs — when trees dominate). Let T(z) be the gener- 
ating function of rooted trees, given by T(z) — ze TtyZ \ Write the generating function of connected 
graphs in the class A as 

T(z) 2 

C(z)=T(z)--f- + D{z). 

The generating function of graphs of A is A{z) = e c ^ z \ As n — > oo, 

e™ 

Cn ~ "! /-— i To and a n ^A(l/e)c n . 

In particular, the probability that Q n is connected tends to l/A(l/e) = e^ 1 ^ 2 ^ 01 -- 1 ^^ as n — > oo. 

Proof. As in Example [2j we use singularity analysis |19l Chap. VI]. By assumption, D(z) has 
radius of convergence larger than 1/e, and the singular behaviour of C{z) is that of unrooted 




FIGURE 3. A random graph of size n = 541 avoiding the diamond, the bowtie 
and the 20-spoon. 
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trees. More precisely, it follows from (|6| that, as z approaches 1/e, 

C{z) = 1/2 + D(l/e) - (1 - ze)(l + D'(l/e)/e) + ^(1 - zef' 2 + 0((1 - zef ), 

o 

this expansion being valid in a A-domain. This gives the estimate of c„ via singularity analysis. 
For the series A, we find 

A{z) = e V2+c(i/ e ) ^1 - (1 - ze ){\ + D'(l/e)/e) + ^(l - ze) 3 / 2 + 0((1 - ze) 2 )^J , 

and the estimate of a n follows. a 

Proposition 7 (Number of components — when trees dominate). The mean of N n 
satisfies: 

E(N n ) ~ 1 + C(l/e) 

and iV n — 1 converges in law to a Poisson distribution of parameter C(l/e) (see ^f§). 

Proof. We can start from (JTJ and apply singularity analysis. Or we can apply a ready-to-use 
result of Bell et al. [H Thm. 2], which uses the facts (proved in Proposition [6]) that the sequences 
c„_i/c„ and c n /a n converge. a 

Proposition 8 (Size of components — when trees dominate). The random variable n—S n 
converges to a discrete limit distribution X given by 

p <* = fc > = i(^)^ 

where ak and A{z) are given in Proposition^ The same holds for n — L n . 

Proof. The two ingredients used in the proof of Proposition [5] to establish the limit law of n — S n 
(namely, smoothness of a n /n\ and convergence of c n /a n ), still hold here (see Proposition [6]) . 
Lemma [4] gives then the law of n — L n . a 

We now present a collection of classes dominated by trees. 

Proposition 9. Let k > 1. Let A be a decomposable class of graphs that contains all trees, and 
such that all graphs of A avoid the k-spoon (shown in Figure [l]). Then A is dominated by trees, 
and the results of Propositions^ [7] and [8] hold. 

Proof. Clearly, it suffices to prove this proposition when A is exactly the class of graphs avoiding 
the fc-spoon, which we henceforth assume. 

We partition the set C of connected graphs of A into three subsets: the set Co of trees, counted 
by Co = T — T 2 /2 with T = T(z), the set C\ of unicyclic graphs (counted by Ci), and finally the 
set C2 containing graphs with at least two cycles (counted by C2). Hence C = T — T 2 /2 + Ci +C2. 
We will prove that C\ has radius larger than 1/e, and that Ci is entire. 

Let Tfe(z) denote the generating function of rooted trees of height less than k (where by 
convention the tree consisting of a single vertex has height 0). That is, T\(z) = z and for k > 1, 

T k+1 (z) = ze T ^. 

Note that Tk(z) is entire. A unicyclic graph belongs to C if and only if all trees attached to its 
unique cycle have height less than k. The generating function of cycles is given by: 



n>3 



Hence, the basic rules of the symbolic method of [19, Chap. II] give: 

CM = Cyc(T fe ) = \ (log T -L R - T k{ z) - ™f) . (9) 
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Recall from Example [2] that T(z) equals 1 at its unique dominant singularity 1/e. Also, T^{z) < 
T(z) for all z € [0, 1/e] since Tk counts fewer trees than T. In particular, Tit (1/e) < 1 and C\(z) 
has radius of convergence larger than 1/e. 

We now want to prove that C 2 is entire. The (2) -core of a connected graph G is the (possibly 
empty) unique maximal subgraph of minimal degree 2. It can be obtained from G by deleting 
recursively all vertices of degree or 1 (or, in a non-recursive fashion, all dangling trees of G). 
By extension, we call core any connected graph of minimal degree 2. Let C 2 denote the set of 
cores having several cycles and avoiding the fc-spoon, and G2 the associated generating function. 
The inequality 

C 2 (z) < C 2 (T k (z)) 

holds, coefficient by coefficient, because the core of a graph of C 2 has several cycles and avoids 
the fc-spoon. Since T k is entire, it suffices to prove that C 2 is entire. It follows from [6, Thm. 3.1] 
that it suffices to prove that no graph G of C 2 contains a path of length 3fc — 1. So let P = 
(t>o, t>i, • • • , vg) be a path of maximal length in G, and assume that i > 3fc — 1. We will prove 
that G contains the fc-spoon as a minor. Since P is maximal and G is a core, there exist vt and 
Vj, with i > 2 and j < £ — 2, such that the edges {vq, Vi] and {vj, vg} belong to G. 

Hi — I or j = 0, let P be the cycle of G formed of P and the edge {vo, vg}. Let Q be another 
cycle of G. If Q contains at most one vertex of P (Figure |4ja) ) , we find an ^-spoon by deleting 
one edge of P, contracting Q into a 3-cycle and one of the paths joining P to Q into a point. If 
Q contains at least two vertices v a and v\, of P. with a < b (Figure [4jb) ) , we may assume that Q 
consists of the edges {v a , u Q +i}, . . . , Vb} and of a path Q that only meets P at v a and Vb- 

Let R denote the cycle formed of the path Q and the path (vb, Vfc+i, • ■ ■ , ve., Vq, . . . , v a }. Then 
we obtain a p-spoon, with p > [3fc/2] — 1 > fc, by contracting the shortest of the cycles Q and 
R into a 3-cycle and deleting an edge ending at v a from the other. 

Assume now that i < £ and j > 0. Suppose first that i < j (Figure |4jc)). By symmetry, 
we may assume that the cycle P\ = (vq, . . . , v) is shorter than (or equal in length to) the cycle 
P 2 = {vj, ■ ■ ■ , Ve). In particular, i < £/2. Contract Pi into a 3-cycle, and remove the edge {vj, ve} 
from P 2 : this gives a p-spoon with p = £ — i > \£/2] > k. Assume now that j < i (Figure md)). 
Consider the three following paths joining Vi and Vf. (Uj, Vi—\, . . . , Vj), (Vi, Vq, Vi, . . . , Vj) and 
(vi,Vi+i, ...,«£, Vj). Since the sum of the lengths of these paths is I + 2 > 3fc + 1, one of them, 
say (vi,Vo,Vi, . . . ,Vj), has length at least fc + 1. That is, j > fc. Delete from this path the 
edge {vi,vo}, and contract the cycle formed by the other two paths into a 3-cycle: this gives a 
j-spoon, with j > k. ■ 




Figure 4. A core having several cycles and avoiding the fc-spoon cannot contain 
a path of length 3fc — 1. 
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The simplest non-trivial class of graphs satisfying the conditions of Proposition [9] consists of 
graphs avoiding the 1-spoon. By specializing to k = 1 the proof of that proposition, we find 
C\ = Cyc(z) and C2 = (since no core having several cycles avoids the 2-path). Hence 

cw =T (= ,-iai + ifi„ g . 1 



2 V 1 - z 

More generally consider the class A^ of graphs avoiding the fc-spoon, but also the diamond 
and the bowtie (both shown in Figure [T| : excluding the latter two graphs means that no graph 
of C can have several cycles, so that C2 = 0. Hence the proof of Proposition [9] immediately gives 
the following result. 

Proposition 10 (No diamond, bowtie or fc-spoon). Let k > 1. Let T(z) be the generating 
function of rooted trees, given by T(z) = ze T ( z \ and let T k (z) be the generating function of 
rooted trees of height less than k. That is, T\(z) = z and for k > 1, 

T k+1 (z) = ze T ^ . 

Let A^ be the class of graphs avoiding the diamond, the bowtie and the k-spoon. The gener- 
ating function of connected graphs of A^ is 

C( k \z)=T(z)-^^+D ik \z) 

where 

T k (zf 



-T k {z)- 



T k (z) ^ ' 2 

The class A^ is dominated by trees, and the results of Propositions^^and^hold. In particular, 
the probability that a random graph of A^ is connected tends to e~ ci ) ( 1 / e ^ as n — > 00. Since 
T k (l/e) tends to T(l/e) = 1 as k increases, this limit probability tends to 0. 

A random graph of ^4n fc ^ is shown in FigureJ^Jfor k — 20 and n — 541. We have also determined 
the generating function of graphs that avoid the 2-spoon. 

Proposition 11 (No 2-spoon). Let T(z) be the generating function of rooted trees, given by 
T(z) = ze T ^ z \ The generating function of connected graphs avoiding the 2-spoon is 

T(z) 2 

C(z)=T(z)-^+D(z) 



where 



The class of graphs avoiding the 2-spoon is dominated by trees, and the results of Propositions^^ 
and [8] apply. 

Proof. We first follow the proof of Proposition [9J we write C = T — T 2 /2 + C\ + C2, where C\ 
is given by ^ with T k = T2 = ze z , and C2 counts connected graphs having several cycles and 
avoiding the 2-spoon. Note that C\ is the first term in the above expression of D(z). Let us now 
focus on C2. 

In Section [10] below, we study the class of graphs that avoid the bowtie, and in particular 



describe the cores of this class (Proposition 31). Since the bowtie contains the 2-spoon, graphs 



that avoid the 2-spoon avoid the bowtie as well. Hence we will first determine which cores of 
Proposition [31] have several cycles and avoid the 2-spoon, and then check which of their vertices 
can be replaced by a small tree (that is, a tree of height 1) without creating a 2-spoon. 

Clearly, the cores of Proposition [3~l"|that have several cycles are those of Figures [T6| [T7| and [18} 
Among the cores of Figure [T6| only K4 avoids the 2-spoon. Moreover, none of its vertices can be 



replaced by a non-trivial tree. This gives the term z 4 /4! in D(z). Among the cores of Figures 17 



andn8] only the ones drawn on the left-hand sides avoid the 2-spoon. In these cores, only the 
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two vertices of degree at least 3 can be replaced by a small tree. The resulting graphs are shown 
in Figure [5] and give together the contribution 

\{ze z f (e z -l-z)+ l -{ze*Y (e* - 1 - z - 

(again an application of the symbolic method of |19l Chap. II]). The proposition follows. a 




Figure 5. Graphs with several cycles avoiding the 2-spoon. 



5. Excluding the diamond and the bowtie: a logarithmic singularity 

Let A be the class of graphs avoiding the diamond and the bowtie (both shown in Figure [lj . 
These are the graphs whose components have at most one cycle. They were studied a long time 
ago by Wright |34| . and the following result has now become a routine exercise. 

Proposition 12 (The number of graphs avoiding a diamond and a bowtie). Let T(z) 
be the generating function of rooted trees, defined by T(z) — ze T ^ . The generating function of 
connected graphs of A is 




Figure 6. A random graph of size n = 859 avoiding the diamond and the bowtie. 
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The generating function of graphs of A is A(z) = e c ^ z \ As n — » oo, 

,e n , 1 e r 



n 



An " ' (2e) 1 /4r(l/4) n 3 / 4 ' 

In particular, the probability that Q n is connected tends to ai speed n -1 / 4 as n — > oo. 



(10) 



Proof. The expression of C(z) is obtained by taking the limit k — > oo in Proposition 10 

We now estimate c n and a n via singularity analysis |19l Sect. VI. 4]. Recall from Example [2] 
that T(z) has a unique dominant singularity, at z = 1/e, with a singular expansion ([6| valid in 
a A-domain. Thus 1/e is also the unique dominant singularity of C(z) and A{z), and we have, 
in a A-domain, 

The asymptotic estimates of c„ and a„ follow. B 

Proposition 13 (Number of components — no bowtie nor diamond). The mean and 
variance of N n satisfy: 

E(N n ) ~logn/4, V(iV n )~log7i/4, 

and the random variable N ". io & n /A convenes in law to a standard normal distribution. 

V lo s«/ 4 

Proof. Using (|6j, the estimate ( |11[ ) can be refined into 

C(«) = ibg( I -^) +A + 0(yi^e), (12) 

where A is a constant, and the proposition is a direct application of |19[ Prop. IX. 14, p. 670]. ■ 

The number of connected components is about 1/4 log n. However, the size of the root com- 
ponent is found to be of order n. More precisely, we have the following result. 

Proposition 14 (Size of the root component — no bowtie nor diamond). The normal- 
ized variable S n /n converges in distribution to a beta law, of density (1 — x)~ 3 / 4 /4 on [0, 1]. In 
fact, a local limit law holds: for x € (0, 1) and k = \_xn\ , 

nF(S n = k) ^ ±(1 - x)- 3 /\ 

The convergence of moments holds as well: for i > 0, 

r(5/4)i! < 



r(i + 5/4) 

Proof. Recall that the existence of a local limit law implies the existence of a global one [HI 
Thm. 3.3]. Thus it suffices to prove the local limit law. But this is easy, starting from the 
rightmost expression in and using (10 I. 

For the moments, let us start from (IJT Our first task is to obtain an estimate of C^ t+1 \z) 
near 1/e. Combining (pi and QH Thm. VI. 8, p. 419] gives, for i > 1, 



C {l+l) {z) 



l / „ \ «+i 



4 V 1 - ze 



We multiply this by the estimate (11) of A(z). apply singularity analysis, and finally use (10 1 
to obtain the asymptotic behaviour of the i^ moment of S n . Since these moments characterize 
the above beta distribution, we conclude [191 Thm. C.2] that S n /n converges in law to this 
distribution. B 

We conclude with the law of the size of the largest component, which we derive from gen- 
eral results dealing with components of logarithmic structures [2\. The following proposition is 
illustrated by Figure [7J 
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FIGURE 7. The distribution function of L n for n = 100. The first plot shows 
the change of regime at n/2, the second the change at n/3, the third the change 
at n/4. 



Proposition 15 (Size of the largest component — no bowtie nor diamond). The nor- 
malized variable L n /n converges in law to the first component of a Poisson-Dirichlet distribution 
of parameter 1/4: for x G (0, 1), 

P(L„ < xn) -> p{l/x), 

where p : R + — > [0, 1] is the unique continuous function such that p(x) = 1 for x G [0, 1] and for 
x > I, 

x x ' A p'{x) + \{x- l)~ 3/4 p(x - 1) = 0. 

A local limit law also holds: for x G (0, 1) and 1/x ^ N, 

(l-x)- 3 / 4 fl-x 

n¥(L n = P 1 

Ax 

Proof. A decomposable class of graphs A is an assembly in the sense of |2] Sec. 2.2] . In particular, 
it satisfies the conditioning relation [5J Eq. (3.1)]. When A is the class of graphs avoiding 



the diamond and the bowtie, ( 10 1 shows that this assembly is logarithmic in the sense of [H 
Eq. (2.15)]: indeed, [2] Eq. (2.16)] holds with irii — Ci, y — e and 9 = 1/4. Our random variable 
L n coincides with the random variable of We then apply Theorem 6.12 and Theorem 6.8 
of this gives the convergence in law of L n and the local limit law. The distribution function 
of the limit law is given by (2j Eq. (5.29)], and the differential equation satisfied by p follows 
from H Eq. (4.23)]. ' a 



Remark. If we push further the singular expansion (12 1 of C(z), we find a term in \/l — ze, 
but its influence is never felt in our asymptotics results. We would obtain the same results (with 
possibly different constants) for any C(z) having a purely logarithmic singularity. 

6. Hayman admissibility and extensions 

Our next examples (Sections [7] to 10 1 deal with examples where C(z) diverges at p with 
an algebraic singularity. This results in A(z) diverging rapidly at p, and leads us to estimate 
a n using a saddle point method, or rather, a black box that applies to Hayman- admissible (or 
H-admissible) functions. Let us first recall what this black box does [19] Thm. VIII. 4, p. 565]. 
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Theorem 16. Let A{z) be a power series with real coefficients and radius p € (0, oo]. Assume 
A{r) is positive for r G (i?, p), for some R £ (0, p). Let 



a\r) = r and byr) = r- 



,A"(r) 
Mr) 



A'(r) 
A(r) 



A{r) v ' A{r) 

Assume that the following three properties hold: 
Hi [Capture condition] 

lim a(r) = lim b(r) = +oo. 

r— yp r^p 

H 2 [Locality condition] For some function 0o(r) defined over (R,p), and satisfying < 
0o(r) < 7r, one has, as r — > p. 



sup 

|e|<9o(r) 



A(re 10 ) i0g( r ) + 6 2 b(r)/2 

A(r) 



1 



0. 



H3 [Decay condition] As 



sup 



A{re ie ) 



A{r) 



Then the rfi 1 coefficient of A(z) satisfies, as n — > oo, 



[z n ]A(z) 



MO 



where C = C„ is the unique solution in (R, p) of the saddle point equation £A'(£) = nA(Q). 

Conditions H 2 and H 3 are usually stated in terms of uniform equivalence as r — > p, but we 
find the above formulation more explicit. 

The set of H-admissible series has several useful closure properties JXSJ Thm. VIII. 5, p. 568]. 
Here is one that we were not able to find in the literature. 

Theorem 17. Let A{z) = F(z)G(z) where F(z) and G(z) are power series with real coefficients 
and radii of convergence < pf < Pg < 00. Assume that F(z) has non-negative coefficients and 
is Hayman-admissible, and that G{pp) > 0. Then A(z) is Hayman-admissible. 

Proof. Let us first prove that the radius of convergence p of A(z) is pp. Clearly, p > pp. 
Now, suppose p > pp. Then A(z) is analytic at pp. Together with G(pp) > this implies 
that F(z) = A(z)/G(z) has an analytic continuation at pp, which is impossible by Pringsheim's 
Theorem (since F(z) has non-negative coefficients). Note also that A(r) is positive on an interval 



of the form [R, p) (by continuity of G). Let us now check the three conditions of Theorem 16 
We have 

a{r) = a F (r) + a G (r), b(r) = b F (r) + b G (r), 

with 



F'(r) 

ap(r) = r and bp(r) 



F'(r) , 2 F"{r) 



F(r) 



F(r) 



F'(r) 
F(r) 



F(r) 

and similarly for oq and bo. 

Hi The capture condition holds for A since it holds for F, given that G(p) > and pc > p. 

H 2 Choose O (r) = 0^ (r) where 0^ (r) is a function for which F(z) satisfies H 2 and H3. We 
have 

A(re l6 ) la ( r )e+e 2 b(r)/2 F(re l6 ) taF ( r )e+e 2 b F (r)/2 G(re l6 ) i aG (r)e+0 2 b 2 a (r)/2 n o\ 
Mr) F(r) G(r) ' ( ' 

By assumption, F satisfies the locality condition: hence 



F(r) 



= 1 + M(r,e) 



(14) 
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where 



sup |M(r,0)| 

\9\<e (r) 



(15) 



as r — > p. For r € [R, p) and \9\ < 0o(r), let us expand log G(re 10 ) in powers of 0: 



log G(re* e ) = log G(r) + i6a G (r) - -b G (r) + 6 6 S(r, 9) 

where S (r, 9) is bounded uniformly in a neighborhood of (p, 0) . We can assume that 9o (r) — > 
as r -> p (see [23l Eq. (12.1)]). Thus 



G ( re ' g ) . p -*aG(r)e+e 2 6 G (r)/2 _ 9 3 S{r,6) _ , N ( 

G(r) ~ i+^r, 



where 



sup \N(r,6)\->0 
\e\<e (r) 



(16) 
(17) 



as r — > p. Putting together Eqs. (13) to (|17l, we obtain that A(z) satisfies H2. 
H 3 We have: 

A (re 10 ) F(re* e )G(re w ) / —- — — 

' - F(r)G(r) V MO + MO 



A(r) 



< 



for r close to p, 



F(r) v w G(r) 

because bp{r) — > 00 as r — ¥ p while bc(r) is bounded around p. Also, since G has radius larger 
than p and G(p) > 0, the term G(re l0 )/G(r) is uniformly bounded in a neighborhood of the 
circle of radius p. Since by assumption, F(z) satisfies H3, this shows that A(z) satisfies it as 
well. a 

We will also need a uniform version of Hayman-admissibility for series of the form e uC ( z K 

Theorem 18. Let C(z) be a power series with non-negative coefficients and radius of conver- 
gence p. Assume A(z) = e c ^ has radius p and is Hayman- admissible. Define 



b{r) = rC'{r) + r 2 C"{r) and V(r) = G(r) 



(rG'(r)) 2 



rC'(r) +r 2 C"(r)' 

Assume that, as r — > p, 

V(r) -> +00, 

-°^L _> 

K ( r )3/a ^ u ' 
^r) 1 /^^ = 0(1). 

Then A(z,u) := e"* 7 ^ safe/^es Conditions (l)-(6), (8) and (9) 0/ [171 Def. 1]. // iV n is a 
sequence of random variables such that 



(18) 
(19) 

(20) 



F(N n = k) 



\z n \e c ^ ' 



</ien </ie mean and variance of N n satisfy: 

E(N n ) ~ G(C„), V(JV„) ~ nC«), (21) 

where ( n = £ is the unique solution in (0, p) of the saddle point equation £C"(£) = n. Moreover, 
the normalized version of N n converges in law to a standard normal distribution: 

N n -E(N n ) 



■A/-(0,1). 



16 



M. BOUSQUET-MELOU AND K. WELLER 



Proof. With the notation of [T71 Def. 1], we have 

a(r, u) — c(r, u) = ruC'(r) — ua(r), b(r, u) = ruC'(r) + r 2 uC"(r) — ub(r), 

K (22) 



a(r,u) — b(r,u) — uC(r), 



e(r) 



for a fixed constant K. Condition (1) of |17l Def. 1] holds for R = p, Rq = and any £. 

idition (4) holds 
p, Condition (5) 



Conditions (8) and (9) hold by the above assumptions (18 1 and (19 1. Condition (4) holds 
because A is Hayman admissible. Since b(r,u)/b(r, 1) = u and e(r) — > as r 



holds. Since a(r,u) = a(r, 1) + c(r, l)(u — 1), Condition (6) obviously holds. 

We are thus left with Conditions (2) and (3) , which are uniform versions (in u) of the locality 
and decay conditions H2 and H3 defining Hayman admissibility. We begin with Condition (2). 
Since A(z) is H-admissible, let 9o(r) be a function for which H2 (and H3) holds: 

A k^l e -iMr)+eH{r)/2 = 1 + M , e) 

A(r) 



where 



M{r) := sup \M{r,( 
\e\<e ( r ) 



tends to as r — > p. Then, for u € [1 — s(r), 1 + £(r)], 

A(re* 9 ,u) e _ iea{rjU)+eH(riU)/2 = ^ ^ ( C ( reid) _ c{r) _ iQa{r) + e%( r )/2)) 
Ayr, u) 

= (1 + M(r,e)) u 

where we have taken the principal determination of log to define (1 + M(r, 6)) u = exp(ulog(l 
M{r,6)) (because M(r,9) is close to 0). Thus 



sup 

|8|<fl (r) 
|«-l|<s(r) 



A(re lS , u) c -ie a (r,u)+6 2 b(r,u)/2 _ ^ 

A(r, u) 



sup |(1 + M(r,6)Y 



II 



e|<e ('') 

«-ll<e(r) 



< 



(1 + e(r))M(r) + O ((1 + e{r))M{rf) , 

and this upper bound tends to as r — > p. This proves Condition (2) with S(r) = 9o(r). 

We finally address Condition (3). Since A(z) satisfies the decay condition H3. the quantity 



N(r, 9) 



A{re l9 ) 
A{r)^bjr) 



satisfies 



sup \N(r,9)\^Q 



(23) 



as r — > p. We have, for u G [1 — s(r). 1 + e(r)] 

A(re l9 ,u) 



A(r, u)y/b(r,u) 



= \N(r,9)\ u 



■e(r) 



< \N(r,9)\^^) 

v 1 - £ ( r ) 



and this tends to uniformly for |0| G [9 {r),n] thanks to (|23j, (|22|, (|18) and ([20 1 . 

As explained in [17] just below Theorem 2, these eight conditions give the estimates (21) of 
E(iV„) and V(N n ) and imply the existence of a gaussian limit law. a 



We finish this section with a simple but useful result of products of series Thm. 2\. 



ASYMPTOTIC PROPERTIES OF SOME MINOR-CLOSED CLASSES OF GRAPHS 



17 



Proposition 19. Let F(z) = ^2 f n z n and G(z) = "Yli n g n z n be power series with radii of 
convergence < pf < Pg < o° ; respectively. Suppose that G(pf) and the sequence (f n )n>0 
is smooth (meaning that f n -i/ fn approaches pp as n — > oo). Then [z n ]F(z)G{z) ~ G{pF)f n - 



7. Graphs with bounded components: C(z) is a polynomial 

Fix k > 1. Let A be the class of graphs in which all components have size at most k. This 
class is minor-closed, with excluded minors all connected graphs of size k + 1. We begin with the 
enumeration of these graphs. The following proposition is a bit more precise than the standard 
result on exponentials of polynomials [191 Cor. VIII. 2, p. 568], since it makes the behaviour of 
b(Cn) explicit. 

Proposition 20 (The number of graphs with small components). The generating func- 
tion of connected graphs of A is the following polynomial: 

i=0 

where Ci denotes the number of connected graphs on i vertices. The generating function of graphs 
of A is A(z) = e c< - z \ As n -> oo, 



^(0 



(24) 



where £ = („ is defined by QC (£) = n and satisfies 

Q = an 1 / k + P + Oin- 1 ^) (25) 

with 

and = J±^±. (26) 



Cfc / fcc fc 
The probability that Q n is connected is of course zero as soon as n > k 



Proof. The series A(z) is H-admissible ([19, Thm. VIII. 5, p. 568]) and Theorem 16 applies. 
The saddle point equation £C"(£) = n is an irreducible polynomial in Q and n, of degree k in 




Figure 8. A random graph of size n = 1171 with component size at most 3. 
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C- Consider 1/n as a small parameter x. By |33l Prop. 6.1.6], the saddle point ( admits an 
expansion of the form 

C=J2 a * n ~ l/k > (27) 

i>i 

for some integer io and complex coefficients on. Using Newton's polygon method, one easily 



finds io = — 1 and the values ( 26 1 of the first two coefficients 



Since b(r) = rC'(r) + r C"(r) has leading term fccfcr /(k — 1)!, the first order expansion of 
6(C) reads 

6(C) = kn + Oin^-^^), 
and the asymptotic behaviour of a n follows. a 

Again, the following proposition is more precise that the statement found, for instance, in [T2l 
Thm. I], because our estimates of E(7V„) and V(N n ) are explicit. Note in particular that E(jV„) ~ 
n/k suggests that most components have maximal size fc. 

Proposition 21 (Number of components — Graphs with small components). The 

mean and variance of N n satisfy: 

E(N n ) ~ Y(N n ) ~ ^a*-i„(fc-i)/*, 



where a is given by (26 1, and the random variable ^==i=p converges in law to a standard 



normal distribution. 



Proof. We apply [121 Thm. I] (we can also apply Theorem 18 if fc > 3). Still denoting the saddle 
point by C = Cn, we just have to find estimates of 

M» = C(C) and a 2 n = C(0 <^>>- 



CC(C) + C 2 C"(C)' 



Given (|25|, we obtain 

= /,■ /,■! 

(fc - 2) 



C 2 C"(C) = (k - l)n - 7i^ ? a fc - 1 n( fc - 1 )/ fe + 0(n^/ k ), 



and finally 



n + C 2 C"(C) ~ fc.fc! 1 j ' 



Since there are approximately n/fc components, one expects the size 5 n of the root component 
to be fc. This is indeed the case, as illustrated in Figure [8] 

Proposition 22 (Size of the components — Graphs with small components). The 

distribution of S n converges to a Dirac law at fc: 



HSU = ]) 



1 ifj = k, 
otherwise. 



The same holds for the size L n of the largest component. 



Proof. We combine the second formulation in d2| with the estimate (24 1 of a n . This gives 



n(j-l)\ A(( n ) Qz] 
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Clearly, it suffices to prove that this probability tends to if j < k. So let us assume j < k. 
Since is increasing with n, it suffices to prove that 

1 C n 



>>n— 7 



0. 



(28) 



Recall from ( 25 1 and ( 27 ) that admits an expansion of the form 

fc-i 

C„=an 1 / fc +^n- J / fe A + 0(l/n). 



i=0 



This gives, for some constants 7j, 

k 

nlogCn = -r logn + nloga + y2j i n 1 ~ i/k + 0(n~ 1/k ) 

Hence 



i=l 



(n- j)logCn-j 



n- 3 



logn + (n- j) loga + ^7 i n 1 ~ i/fe + 0(1). 



i=l 



This gives 



j -k 



n log („ - (ra - j) log £„_j - log n = — — log n + 0(1), 



and (28 1 follows, since j < k. Since L n > SVi, the behaviour of L n is then clear. a 

8. Forests of paths or caterpillars: a simple pole in C(z) 

Let .A be a decomposable class (for instance defined by excluding connected minors), with 
generating function A(z) = exp(C(z)). Assume that 

Q 



C{z) 



D(z), 



1-z/p 

where D has radius of convergence larger than p. Of course, we assume a > 0. 



(29) 



Proposition 23 (The number of graphs — when C has a simple pole). Assume the 
above conditions hold, and let (3 = D(p). As n —¥ oo, 

a l/4 e «/2+£ 

c n ~n\ap " and a„ ~ n! p n e v™. (30) 

In particular, the probability that Q n is connected tends to at snppd. n 3 / 4 p~ 2 ^*™ 



FIGURE 9. A random forest of paths of size n = 636 (left) and a forest of 
caterpillars of size n — 486 (right). 



20 



M. BOUSQUET-MELOU AND K. WELLER 



Proof. The asymptotic behaviour of c„ follows from |19[ Thm. IV. 10, p. 258]. To obtain the 
asymptotic behaviour of a n , we first write 



A(z) = F(z)G(z) with F(z) = exp 



and G{z) = 



- p D(z) 



(31) 



1 - z/p, 

where G(z) has radius larger than p. To estimate the coefficients of F, we apply the ready-to- 
use results of Macintyre and Wilson [26, Eqs. (10)-(14)], according to which, for a, 7 > and a 
non-negative integer k, 



log 



1 



1 



l-z) (1-z)' 



exp 



a 
1 - z 



l/4 p a/2 /n y/2 /log 71 



a > e 



(32) 



This gives 



Q,l/4 e «/2 

/„ := [z n ]F(z) ~ ^^p-e 2 ^. 



This shows in particular that f n -i/ fn tends to p as n —¥ 00, so that we can apply Proposition 19 
to (31 1 and conclude. a 



Proposition 24 (Number of components — when C has a simple pole). Assume (29 1 
holds. The mean and variance of N n satisfy: 



E(N n 



V(iV n 



1/2, 



and the random variable converges in law to a standard normal distribution. 



Proof. We apply Theorem 18 The H-admissibility of A(z) follows from Theorem 17 using (31 I 
and the H-admissibility of exp(a/(l — z)) (see [19, p. 562]). Conditions ( 18 1— ( 20 1 are then readily 
checked, using 

C(r) ~ — , b(r) ~ i—r^ and V(r) ~ — . 

y ' l-r/p' v ; (1-r/p) 3 v ; 2(1- r/p) 

We thus conclude that the normalized version of N n converges in law to a standard normal 
distribution. For the asymptotic estimates of E(N n ) and Y(N n ), we use (21) with the saddle 
point estimate = p — pWa/n + 0(l/n). a 

Since there are approximately y/n components, one may expect the size S n of the root com- 
ponent to be of the order of y/n. 

Proposition 25 (Size of the root component — when C has a simple pole). The 

normalized variable S n / \JnJa converges in distribution to a Gamma(2) law, of density xe~ x on 
[0,oo). In fact, a local limit law holds: for x > and k — \_xy/n/a\, 

\AVa ¥(S n = k) -> xe~ x . 
The convergence of moments holds as well: for i > 0, 

nsi)~(i+mn/ay/ 2 . 

Proof. For the local (and hence global) limit law, we simply combine ^ with (30 1. For the 
moments, we start from (El, with 



C {t+1) {z) 



a(i + l)\ 



p i+1 (l- z/ P y+ 2 



D^ +1 \z). 



Let us first observe that (30) implies that a n /n\ is smooth. We can thus apply Proposition 19 
to the product D^ l+1 \z)A(z), which gives 

[z n - l - 1 }D (t+1 \z)A(z) D( i+1 \p) a * r><t+i)(n\ . 
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We thus have 



(n-1)! 



E(#) 



a(i + l)! 



y +i (i - z/p)^ 



exp 



1 - z/p 



D{z) 



Now (32 1 gives 



J p l+1 (l-^/p)' i+2 



exp 



1-z/p 



a V4 e a/2 , 

a(z + l) — = q; , - 



■/2+1 



(33) 



(34) 



In particular, this sequence of coefficients is smooth. Hence by Proposition [19] the asymptotic 
behaviour of (33) only differs from ( 34 1 by a factor e' 3 , where j3 = D(p). Combined with (30), 



this gives the limiting i^ 1 moment of S n . Since these moments characterize the above Gamma 
distribution, we can conclude |19| Thm. C.2] that S n / \fnfa converges in law to this distribution. 



We now present two classes for which C(z) has a simple isolated pole: forests of paths, and 
forests of caterpillars (a caterpillar is a tree made of a simple path to which leaves are attached, 
see Figure [T]). In forests of paths, the excluded minors are the triangle K% and the 3-star. The 
fact that N n converges in probability to y/n/2 for this class was stated in [27, p. 587]. For forests 
of caterpillars, the excluded minors are K3 and the tree shown in Table ll| (6 
is also considered in [B]. It is also the class of graphs of pathwidth 1. 



Proposition 26 (Forests of paths or caterpillars) 

of caterpillars are respectively 



C p {z) 



2(1 -z) 



and 



C c {z) 



th line). This class 

The generating functions of paths and 
If 



2(1 



ze^ 



z 
~2 



(35) 



For both series, Condition (29 I is satisfied and Propositions 23 24 and 25 hold. For paths we 
have p = 1, a = 1/2 and f3 := D(p) — 0. For caterpillars, p ~ 0.567 is the only real such that 
pe" = l, 



a 



0.06 and f3 



p (l0 + 3p-4p 2 



2(1+ p) 4(l + p) / 

Proof. The expression of C p (z) is straightforward. One can also write 

C p (z) 



0.59. 



(36) 



1 



1 



2(1 -z) 



which gives D p (l) = 0. Let us now focus on caterpillars. Let us call star a tree in which all 
vertices, except maybe one, have degree 1. By a rooted star we mean a star with a marked vertex 
of maximal degree: hence the root has degree 1 for a star with 1 or 2 vertices, and degree at 
least 2 otherwise. Clearly, there are n rooted stars on n labelled vertices, so that their generating 
function is 

S'(z) = > 7 — , = ze z . 

w ^ (n - 1) 

n>l V ' 

The generating function of (unrooted) stars is 

v 2 Jl 



S(z) 



Z" 



(because all stars have only one rooting, except the star on 2 vertices which has two). Now a 
caterpillar is either a star, or is a (non-oriented) chain of at least two rooted stars, the first and 
last having at least 2 vertices each. This gives 

(S'(z)-z) 2 



C(z) = S{z) + 



2(1- S'(z)Y 



which is equivalent to the right-hand side of (35 
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The series C c {z) is meromorphic on C, with a unique dominant pole at p, and its behaviour 
around this point is easily found using a local expansion of ze z at p: 

C c (z) = -^— +/3 + 0(l-z/p), 
1 - z/p 



with a and /? as in ( 36 1 . ■ 

For forests of paths, we have also obtained the limit law of the size L n of the largest component. 
It is significantly larger than the root component (y^logn instead of y/n). 

Proposition 27 (Size of the largest component — forests of paths). In forests of paths, 
the (normalized) size of the largest component converges in law to a Gumbel distribution: for 
x € K and as n — > oo , 



L„ — y/n/2 log n 
y/n/2 



< x I — > exp 



,-x/2 



Proof. We start from ([5]), where 

k — \Jn/2 (log n + x) 
and the generating function of paths of size less than k is: 



C [k] (z) = 



2 ' 2(1-2) 

Using a saddle point approach for integrals [191 p. 552], we will find an estimate of 



[z \ exp 



(CM(z)) 



1 

2lTT 



exp 



(CW(*)) 



2: 



dz 

n+l 1 



(37) 
(38) 

(39) 



where the integration contour is any circle C r of center and radius r < 1. 

Let us first introduce some notation: we denote C^(z) by K (z), the integrand in (39 1 by F, 
and its logarithm by /: 

K(z) = CW(z), F(z) = eM * + [ Z)) , f(z)=K(z)-(n + l)\ogz. 

We choose the radius r = r n that satisfies the saddle point equation 

F'(r) = 0, or equivalently f'{r) — or rK'(r) = n + l. 

Note that rK'(r) increases from to oo as r grows from to 1, so that the solution of this 
equation is unique and simple to approximate via bootstrapping. We find: 



r = 1 - 



1 



"*/2 bp 



2n 4y/2 



O 



(40) 



Gaussian approximation. Let 0q £ (0, 7r). By expanding the function g : 9 n- f(re zS ) in the 
neighbourhood of 6> = 0, we find, for \6\ < 9q: 



\f{re w ) - f(r) + 6 2 r 2 f(r)/2\ < f sup g^{a) 



() 



|a|<e 



(41) 



with 



,(3) 



(a) 



< K'(r) + 



ire la f (re la ) - 3ir 2 e Ma f" (re la ) - ir A e Ma f" (re la ) \ 
n + l „ T ,„, , „n+l T , „n+l 



- + 3K"(r) + 3^ + K"'{r) + 2™ 



By combining the expression (38 1 of K(z) = C^- k \z) and the saddle point estimate (40), we find 
that K'(r) = (n+ l)/r - n, that K"(r) - 2V2n 3 / 2 and finally that K"'{r) - 12n 2 . This term 
dominates the above bound on |g( 3 )(a)|. Hence, if 

e = 9 Q (n)=o(n-^ 3 ), (42) 
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we find, by taking the exponential of ( |4lj ), 

F(re M ) ~ F(r)e- e2r2 f"^/ 2 , 

uniformly in \6\ < 8 . 

Completion of the Gaussian integral. We split the integral (39 1 into two parts, depending 
on whether \9\ < da or \0\ > da- The first part is 

,re lB dd rF(r) 



e 



F(re w )- 



2tt 



2tt 



-e 2 r*f"(r)/2 e i6 d ^ 



On 



As argued above, r 2 /"(r) ~ K" (r) ~ 2\[2r?l' 1 . Hence, if we choose 6q = 6o(n) such that 
#q7j 3 / 2 — > oo (which is compatible with (42), for instance if 

#o - n- 5 / 7 , (43) 

which we henceforth assume), we obtain 



F{re i6 )- 



3 dd 



F(r) 



2tt 



e ry/f"(r) 



-a 2 /2 e i a /r/^f"(r) da 



F(r) 



F(r) 
V27r/"(r) 

2 1 /42 v ^n 3 /' 



exp 



,-x/a 

"7T 



(44) 



by gOl. 

The second part of the integral can be neglected. The second part of the integral ( 39 1 is 

re i9 d6 



F(re l 



o <\e\<7T 



2ir ' 



and we want to prove that it is dominated by (44). It suffices to prove that for 9o < \9\ < ir, 

/ F(r) \ 



\F{re w )\=o 

Let us denote z = re 10 and za = re . We have 
\F(re ie )\ 



(45) 



F(r) 



= \exp(K(z)-K(r))\ 

< exp(\K{z)\ - K(r)) = exp 

" exp l2^-2(T37) 

< exp 



2(1 -z) 



2(1 -r) 



2|l-z | 2(1 -r) 



exp 



a l/14 



(l + o(l)) , 



given the values (|37|, (|40} and Q of fe, r and O . Since /"(r) ~ 2V^n 3/2 , we conclude that ([45 1 
holds. 



Conclusion. We have now established that the integral (39) is dominated by its first part, and 
is thus equivalent to (44 1. To obtain the limiting distribution function, it remains to divide this 
estimate by a n /n\. The asymptotic behaviour of a n is given by (30 1, with a — 1/2, p — 1 and 
(3 = 0, and this concludes the proof. ■ 
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Figure 10. A random graph of size n — 1034 avoiding the 3-star. 

9. Graphs with maximal degree 2: a simple pole and a logarithm in C{z) 

Let A be the class of graphs of maximal degree 2, or equivalently, the class of graphs avoiding 
the 3-star. The connected components of such graphs are paths or cycles. This class differs from 
those studied in the previous section in that the series C(z) has now, in addition to a simple 
pole, a logarithmic singularity at its radius of convergence p. As we shall see, the logarithm 
changes the asymptotic behaviour of the numbers a n , but the other results remain unaffected. 
The proofs are very similar to those of the previous section. 

Proposition 28 (The number of graphs of maximal degree 2). The number of connected 
graphs (paths or cycles) of size n in the class A is c n = nl/2 + (n — l)!/2 for n > 3 (with 
C\ = C2 = 1) and the associated generating function is 

, z(2-z + z 2 ) 1 1 



4(1 -z) 2 °l-z 
The generating function of graphs of A is 



A(z) = e° [z > = —j= exp 



yr^^'V 4(i - Z ) 

As n — > oo, 

1 



V2n_ 



™ o / 1 II 

2y/enn 1 ' z 

In particular, the probability that Q n is connected tends to at speed 7 j 1 / 2 e -v / 2n as n ^ oo. 

Proof. Again, the exact results are elementary. To obtain the asymptotic behaviour of a n , we 
write 

A{z) = F(z)G(z) with F(z) = exp ^ and G(z) = exp (~ - ? - 

(46) 



and combine Proposition 19 with (32) 



For the number of components, we find the same behaviour as in the case of a simple pole 



(Proposition 24 with a — 1/2). We have also determined the expected number of cyclic compo- 
nents. 

Proposition 29 (Number and nature of components — Graphs of maximal degree 

2). The mean and variance of N n satisfy: 

E(N n ) ~ V(N n ) ~ ^V8, 
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and the random variable 



N„ 



ln/2 



(n/8) 1 / 4 converges in law to a standard normal distribution. 
The expected number of cycles in Q n is equivalent to (logn)/4. 



Proof. We want to apply Theorem 18 To prove that A{z) is Hayman-admissible, we apply 



Theorem 17 to (46). This reduces our task to proving that F(z) is H-admissible, which is 
done along the same lines as [191 Ex. VIII. 7, p. 562] (see also the footnote of |2"3"l p. 92]). 
Conditions (18 1— (201 are readily checked. The asymptotic estimates of E(iV„) and V(7V„) are 
obtained through (21), using the saddle point estimate £ n = 1 — 1/V271+ 0(l/n). 



The bivariate generating function of graphs of A, counted by the size (variable z) and the 
number of cycles (variables v) is 

•*2 



A(z,v 



exp z 



+ v Cyc(z) 



where Cyc(z) is given by 
in Q n is found to be: 



2(1 -z) 

By differentiating with respect to v, the expected number of cycles 
[z n ] Cyc{z)A{z) 



z n )A(z) 



The asymptotic behaviour of [,z n ]A(z) = a n /n\ has been established in Proposition 28 We deter- 
mine an estimate of [z n ] Cyc(z)A(z) is a similar fashion, using a combination of Proposition 19 
and ((32j. We find 

logn 



[z n ] Cyc(z)A(z) 



I2n 



enn 



1/2 



and the result follows. 



The size of the root component is still described by Proposition 25 with a = 1/2. The proof 
is very similar, with now 



C {l+1 \z) = 



i\ 2 + i - z 

2(i- z y+ 2 



i 



10. Excluding the bowtie: a singularity in (1 — z/p) ~ 1 / 2 

We now denote by A the class of graphs avoiding the bowtie (Figure [T]) . 

Proposition 30 (The generating function of graphs avoiding a bowtie). Let T = T(z) 
be the generating function of rooted trees, defined by T(z) — ze T ^ . The generating function of 
connected graphs in the class A is 



C(z) = 



T 2 (l — T + T 2 ) e 1 



^ 0g (l^T 



T(12-54T + 18T 2 -5T 3 - T 4 



l-T 

The generating function of graphs of A is A(z) = e c ( z \ 



24(1 — T) 



(47) 



This is the most delicate enumeration of the paper. The key point is the following character- 
ization of cores (graphs of minimal degree 2) avoiding the bowtie. 

Proposition 31. The cores that avoid the bowtie are: 

• the empty graph, 

• the cycles, 



• K^, with one edge possibly decorated by a chain of vertices of degree 2 (Figure 16), 

• the graphs of Figures |17| anrf |18[ 

We will first establish a number of properties of cores avoiding a bowtie. Recall that a chord 
of a cycle C is an edge, not in C, joining two vertices of C. 
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Figure 11. A random graph of size n = 758 avoiding the bowtie. 



Lemma 32. Let C = (vq, V\, . . . , u n — l) be a cycle in a core G avoiding the bowtie. Every chord 
of C joins vertices that are at distance 2 on C (we say that it is a short chord). Moreover, C 
has at most two chords. If it has two chords, say {vq,v 2 } and {vi,Vi+ 2 }, with 1 < i < n — 1 and 
v n = v , v n+1 = Vi, then v. t = v x or v i+2 = v x . 

Proof. If a chord were not short, contracting it (together with some edges of C) would give a 
bowtie. Figure [12] then proves the second statement, which can be loosely restated as follows: 
the two chords cross and their four endpoints are consecutive on C . ■ 




no no yes 



Figure 12. The relative positions of two (short) chords in a cycle. 

Lemma 33. Let C be a cycle of maximal length in a core G avoiding the bowtie. Let v be an 
external vertex, that is, a vertex not belonging to C. Then v is incident to exactly two edges, 
both ending on C . The endpoints of these edges are at distance 2 on C . 

Proof. Since G is a core, v belongs to a cycle C' . Since C is connected and avoids the bowtie, C' 
shares at least two vertices with C. Thus let Pi and P 2 be two vertex-disjoint paths (taken from 
C") that start from v and end on C without hitting C before. Let v x and v 2 be their respective 
endpoints on C. Then v x and v 2 lie at distance at least 2 on C, otherwise C would not have 
maximal length. The same argument implies that Pi and P 2 have length 1 each. 

Assume now that v has degree at least 3, and let e be a third edge (different from Pi and P 2 ) 
adjacent to v. Again, e must belong to a cycle, sharing at least two vertices with C, and the 
same argument as before shows that e ends on C. But then Figure [13] shows that G contains a 
bowtie. ■ 
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Figure 13. A cycle C with an external vertex v of degree at least 3. 



Lemma 34. Let C be a cycle of maximal length in a core G avoiding the bowtie. If C has two 
chords, it contains all vertices of G . 

Proof. Let e\ and ei be the two chords of C . Lemma |32| describes their relative positions. Let 
v be a vertex not in C. Lemma [33l describes how it is connected to C. Contract one of the two 
edges incident to v to obtain a chord of C. By Lemma [32] this chord must be a copy of ei or 
e2- But then Figure 14 a) shows that G contains a bowtie (delete the two thick edges). a 




(b) 



no 



Figure 14. A cycle C with (a) two chords and an external vertex v (b) one 
chord and an external vertex v. 



Lemma 35. Let C be a cycle of maximal length in a core G avoiding the bowtie. If C has a 
chord e, all external vertices of C are adjacent to the endpoints of e. 



Proof. Let v be external to C. Contract one of the incident edges. This gives a chord 



If 



e' is a copy of e, then we are done. Otherwise, the relative positions of e and e' are described 
by Lemma 32 But then Figure 14 b) shows that C has not maximal length (consider the cycle 



shown with the dotted line) . a 

Lemma 36. Let C be a cycle of maximal length in a core G avoiding the bowtie. If C has 
several external vertices, they are adjacent to the same points of C . 

Proof. Consider two external vertices v\ and i> 2 . Lemma [33] describes how each of them is 
connected to C. Contract an edge incident to v\ and an edge incident to V2- This gives two 
chords of C. Either these two chords are copies of one another, which means that v\ and v 2 are 
adjacent to the same points of C. Or the relative position of these two chords is as described in 



Lemma 32 But then Figure 15 shows that G contains a bowtie (contract e). 




FIGURE 15. A cycle C with two external vertices. 
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Proof of Propositions and [?7J Observe that a graph G avoids the bowtie if and only if its 
core (as defined in the proof of Proposition |9| avoids it. Hence, if C{z) denotes the generating 
function of non-empty cores avoiding the bowtie, we have 

T{zf 



C(z) = T(z) 



C{T{z)). 



(48) 



Using the above lemmas, we can now describe and count non-empty cores avoiding the bowtie. 
We start with cores reduced to a cycle: their contribution to C{z) is given by pj. We now 
consider cores G having several cycles. Let C be a cycle of G of maximal length, chosen so that 
it has the largest possible number of chords. By Lemma[32j this number is 2, 1 or 0. 



If C has two chords, it contains all vertices of G (Lemma 34 1. By Lemma 32 and Figure 12 
(right), either all vertices have degree 3 and G = K4, or G consists of K4 where one edge is 
substituted by a (directed) path (Figure 16 1. 



> 




FIGURE 16. K4 with a decorated edge. 
This gives the generating function 



z 

if 



z a z 
4! - 6 ' — 



(49) 



where, in the second term, we read first the choice of the 4 vertices of degree 3 forming a K4, 
then the choice of one edge of this K 4 , and finally the choice of a (directed) path placed along 
this edge. 

Assume now that C has exactly one chord e. By Lemma [35] all external vertices are adjacent 
to the endpoints of e. To avoid problems with symmetries, we count separately the cores where 



C has length 4, or length > 5 (Figure 17). This gives the generating function 

^2 



-(e z -l-z) + -(e 2 -l) — . 
2 V ' 2 V 'l - z 



(50) 



In the second term, the factor z 2 /(l — z) accounts for the directed chain of vertices of degree 2 
lying on the maximal cycle. 



> 2 





at least two vertices 



at least one vertex 



Figure 17. A maximal cycle with one chord e and external vertices. 



Assume finally that C has no chord. By Lemma |36[ all external vertices are adjacent to the 
same points of C . Again, we treat separately the cases where C has length 4, or length > 5 
(Figure 18 1. This gives the generating function 

^2 



Z -(e*-l-z-z 2 /2)+ Z -(e* 



l-z) 



1 



(51) 



Putting together the contributions (pj, (49 1, (50) and ( pi) gives the value of C(z) (the gen- 
erating function of cores) from which we obtain the series C{z) using (48 1. 1 
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FIGURE 18. A maximal cycle C with no chord and several external vertices. 



We now derive asymptotic results from Proposition [30] 
Proposition 37 (The asymptotic number of graphs avoiding the bowtie). As n — > oo, 

, e - 5/4 e" 



and 



. (e - 5/4)VV 9 / 8 - llc /3 e » , 
a„ - n! ^_ ^ exp -(e - 5/4) 2/3 n 1/3 



()7T 



2/3 



Proof. Let us first recall that the series T(z) has radius 1/e, and can be continued analytically 
on the domain T> :— C \ [1/e, +oo). In fact, T(z) = — W(— z), where W is the (principle branch 
of the) Lambert function |13| . The singular behaviour of T(z) near 1/e is given by (|6|. Moreover, 
the image of I? by T avoids the half-line [1, +oo). 

It thus follows from the expression (47 1 of C(z) that C(z) and A(z) are analytic in the domain 
V. Moreover, we derive from (pi that, as z approaches 1/e in a A-domain, 



C(z) 



e- 5/4 



(52) 



y/2y/l - ze ' 

The above estimate of c„ then follows from singularity analysis. 

We now embark with the estimation of a n . We first prove (see Proposition 40 in Appendix 
A) that A(z) is H-admissible. We then apply Theorem [16] The saddle point equation reads 
CC"(C) = n. Using the singular expansion (JsJ) of T(z), and a similar expansion for T'(z), this 
reads 



e-5/4 



1 



2^2(1 - Ce) 3/2 4 (! ~ Ce) 
This gives the saddle point as 

£_1 (e-5/4) 2 / 3 



1 



1 



(1 - Ce)Va 



+ <3(7i- 4 / 3 ) 



n. 



(53) 



(54) 



We now want to obtain estimates of the values A(C), C™ an d 6(C) occurring in Theorem [l6| 
Refining ( 52 1 into 



™ ^ e-5/4 1 2 

C Z = /o /i — + Z l0g 1 

V2%/1 - ze 4 1 - ze 



we find 



which gives 



19 lie 

~8~ ~ IT 

53 lie 



+ 0(y/l-ze), 
+ 0(n' 1 / 3 ), 



(55) 



A(0 



53 _ lie 
e24 3 



It then follows from ( 54 1 that 

c n - 



,, _. 5/4 ^- 1/b exp((e-5/4) 2 / 3 nV^. 
e" 1 / 6 exp (-n - (e - 5/4) 2 / 3 n 1 ^) . 



(56) 
(57) 
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Finally, 



so that 



b(r) = rC'{r)+r 2 C"{r) 



6(0 



(e- 5/4) 2 / 3 



3V2{e - 5/4) 
8(1 -re) 5 / 2 ' 

,5/3 



(58) 



Putting this estimate together with (56 1 and ( |57[ ), we obtain the estimate of a n /n\ given in the 
proposition. a 

Proposition 38 (Number of components — no bowtie). The mean and variance of N n 
satisfy: 

E(N n ) ~ (e - 5/4) 2 / 3 n x / 3 , V(7V„) ~ ~(e- 5/4) 2 /V/ 3 , 
and the random variable jVn ~ K ( JV " i converqes in law to a standard normal distribution. 



Proof. We want to apply Theorem 18 By Proposition 40 is H-admissible. Conditions ( 18 1- 

V2(e - 5/4) 



(20 1 are readily checked, using 
e- 5/4 



C(z) 



V2yT 



3vT 



ze 



The asymptotic estimates of E(iV n ) and Y(N n ) are obtained through (21 ), using the saddle point 
estimate (54 1. a 

Since there are approximately rt 1 / 3 components, one may expect the size S n of the root 
component to be of the order of rt 2 / 3 . More precisely, we have the following result. 

Proposition 39 (Size of the root component — no bowtie). The normalized variable 
(e — 5/4) 2 / 3 5' Tl /(2n 2 / 3 ) converges in distribution to a Gamma(3/2) law, of density 2y/xe~ x / 'y/n 
on [0, oo). In fact, a local limit law holds: for x > and k = [x t e ^uyt/s J > 



2n 2 / 3 



3 V(Sn = k) 



(e-5/4) 2 / 3 

The convergence of moments holds as well: for i > 0, 



7T 



r(< + 3/2) 



2n 2 / 3 



r(3/2) V(e-5/4) 2 / 2 



Proof. The local (and hence global) limit law follows directly from Proposition 37 usin g p] ) 



For the convergence of the moments, we start from (|4J. We first prove (see Proposition 40 
Appendix A) that C^ +1 - ) (z)A(z) is H-admissible. We then apply Theorem 16 to estimate the 
coefficient of z n in this series (we will replace n ny n — i — 1 later) . Our calculations mimic those 
of Proposition |37| but the saddle point equation now reads 



_ Ai) 

using the estimate (67) of C^'(z), this reads 

e-5/4 



where £ = (A depends on i and rt. Comparing with the original saddle point equation (53 1, and 

;ads 

7 + 4i 



2y/2(l - Ce) 3 / 2 4(1 - Ce) 
This gives the saddle point as 



O 



1 

(1 - Ce) 1 / 2 



c = 



1 (e-5/4) 2 / 3 7 + 4i 



2en 2 / 3 



6en 



+ 0(n- 4 / 3 ). 



(59) 
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We now want to obtain estimates of C^ +1 ^(C)^4(C), C™ an d h(0- We first derive from (55 1 that 
C(C) = ( e -5/4) 2 /V/ 3 + ilog U 

This gives 



29 lie 2i _ 1/3 . 
6-e-5/4 + 24-^T-3 +0( " ' } ' 



A(0 



53 _ lie _ 2 
6 24 3 3 



(e - 5/4)V8 
Moreover, we derive from J671) that 



-^n 1 / 6 exp((e-5/4) 2 / 3 n 1 / 3 ). 



C (l+1) (C) 



(2i + l)! 
2 4 i! 



,l+2i/3 



It then follows from ( 59 1 that 



c 



-7/6-2i/3 



(e-5/4) 2 */ 3 



cxp ( -n - (e - 5/4) 2/3 n 1/3 /2 



(60) 
(61) 
(62) 



Finally, (69) and (58 1 give 



MO ~ 6(0 



(e- 5/4) 2 / 3 



,5/3 



Putting this estimate together with ( 60 1 , ( |61[ ) and (62 1, we obtain 

[z n ]C {l+1) {z)A{z) ~ 

(22 + 1)! (e - 5 /4)l/6 e 19/8-lle/3 e «+ I+ l n l+2i/3 



2*i! 



67T 



j2/3 ( e _5/4)2i/3 



cxp -(e-5/4) 



2/3„l/3 



We finally replace n by n — i — 1 (the only effect is to replace e ™ +l+1 by e n ), and divide by the 
estimate of na n /n\ given in Proposition 37 this gives the estimate of the ft 1 moment of S n as 



stated in the proposition, and concludes the proof. a 

11. Final comments and further questions 

11.1. Random generation 

For each of the classes A that we have studied, we have designed an associated Boltzmann 
sampler, which generates a graph G of A with probability 

where x > is a fixed parameter such that A(x) converges. We refer to |18l Sec. 4] for general 
principles on the construction of exponential Boltzmann samplers, and only describe how we 
have addressed certain specific difficulties. Most of them are related to the fact that our graphs 
are unrooted. 

Trees and forests. Designing a Boltzmann sampler for rooted trees is a basic exercise after 
reading (THj. Note that the calculation of T(x) can be avoided by feeding the sampler directly 
with the parameter t = T(x), taken in (0, 1]. To sample unrooted trees, a first solution is to 
sample a rooted tree G and keep it with probability 1/|G|. However, this sometimes generates 
large rooted trees that are rejected with high probability. A much better solution is presented 
in [T4"l Sec. 2.2.1]. In order to obtain an unrooted tree distributed according to (63 1, one calls 
the sampler of rooted trees with a random parameter t. The density of t must be chosen to be 
(1 - t)/C{x) on [0,T(x)], with C{x) = T(x) - T(x) 2 /2. To sample t according to this density, 
we set t = 1 — Wl — 2uC(x), where u is uniform in [0, 1]. Again, we actually avoid computing 
the series C(x) by feeding directly our sampler with the value T(x) g (0, 1]. We use this trick 
for all classes that involve the series T(x). 

To obtain large forests (Figure [2]), we actually sample forests with a distinguished vertex; that 
is, a rooted tree plus a forest [HI Sec. 6.3]. 
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Paths, cycles and stars. The sequence operator of |18l Sec. 4] produces directed paths, while 
we need undirected paths. Let u be uniform on [0, 1]. Our generator generates the one- vertex 



path if u < x/C p (x), where C p (x) is given by (351, and otherwise generates a path of length 



2 + Gcom(a;). An alternative is to generate a directed path, and reject it with probability 1/2 if 
its size is at least 2. 

Although the cycle operator of |18l Sec. 4] generates oriented cycles, this does not create a 
similar problem for our non-oriented cycles: indeed, a cycle of length at least 3 has exactly two 
possible orientations. 

Designing a Boltzmann sampler TRS for rooted stars is elementary. For unrooted stars, we 
simply call TRS, but reject the star with probability 1/2 if it has size 2 (because the only star 
with two roofings has size 2). 

Graphs avoiding the bowtie. This is the most complex of our algorithms, because the 



generation of connected graphs involves 7 different cases (see the proof of Proposition 30 ) . There 



is otherwise no particular difficulty. We specialize this algorithm to the generation of graphs 



avoiding the 2-spoon (Proposition 11 1. However, the probability to obtain a forest is about 0.95, 
and thus there is no point in drawing a random graph of this class. 

The graphs shown in the paper have been drawn with the graphviz software. 

11.2. The nature of the singularity of C(z) 

This is clearly a crucial point, as the probability that Q n is connected and the quantities 
N n and S n seem to be directly correlated to it (see the summary of our results in Table [l). 
This raises the following question: is it possible to describe an explicit correlation between the 
properties of the excluded minors and the nature of the singularity of C(z)l For instance, it is 
known that C(p) is finite when all excluded minors are 2-connected, but Section [4] shows that 
this happens as well with some non-2-connected excluded minors. Which excluded minors give 
rise to a simple pole in C(z) (as in caterpillars)? or to a logarithmic singularity (as for graphs 
with no bowtie nor diamond), or to a singularity in (1 — z/p) -1 / 2 (as for graphs with no bowtie)? 

Some classes for which C{z) has a unique dominant pole of high order are described in the 
next subsection. 

11.3. More examples and predictions 

Our examples, as well as a quick analysis, lead us to predict the following results when C{z) 
has a unique dominant singularity and a singular behaviour of the form (1 — z/ p)~ a , with a > 0: 

• the mean and variance of the number N n of components scale like n a ^ 1+a \ and N n 
admits a gaussian limit law after normalization, 

• the mean of S n scales like n 1 ^ 1+a \ and S n , normalized by its expectation, converges to 
a Gamma distribution of parameter a + 1 . 

The second point is developed in [11 . To confirm these predictions one could study the following 
classes, which yield series C(z) with a high order dominant pole. Fix k > 2, and consider the 
class of forests of degree at most k, in which each component has at most one vertex of 
degree > 3. This means that the components are stars with long rays and "centers" of degree at 
most k. It is not hard to see that 



k 



2(1 -z) ^i\(i- z y 

so that Ck has a pole of order k (for k > 3). The case k — 2 corresponds to forests of paths 
(Section [8]). The limit case k = oo (forests of stars with long rays) looks interesting, with a very 
fast divergence of C at 1 : 

~2 



C* (oo) (z) = zexp 



2(1 -zf 



We do not dare any prediction here. 
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11.4. Other parameters 

We have focussed in this paper on certain parameters that are well understood when all 
excluded minors are 2-connected. But other parameters — number of edges, size of the largest 2- 
connected component, distribution of vertex degrees — have been investigated in other contexts, 
that sometimes intersect the study of minor-closed classes [H 13 [HI [16j [21] . When specialized to 
the theory of minor-closed classes, these papers generally assume that all excluded minors are 
2-connected, sometimes even 3-connected. 

Clearly, it would not be hard to keep track of the number of edges in our enumerative results. 
Presumably, keeping track of the number of vertices of degree d for any (fixed) d would not be 
too difficult either. This may the topic of future work. The size of the largest component clearly 
needs a further investigation as well. 
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Appendix A: Hayman-admissibility for bowties 



Proposition 40. Let C(z) and A{z) be the series given in Proposition 30 Then the series A{z) 
and C( I+1 \z)A(z) are H-admissible, for any i > 1. 

Proof. We begin with the series A(z). Recall the analytic properties of T(z), listed at the 
beginning of the proof of Proposition |37[ The capture condition Hi is readily checked. In fact, 

e - 5/4 



a(r) = rC'{r) 



2^(1 - re) 3 / 2 



and b(r) = rC'(r) + r 2 C"{r) ~ ^ -^g . (64) 



Let us now prove H 2 . By Taylor's formula applied to the function / : 9 h-> C(re l9 ), we have, 
for r € (0, 1/e) and 9 € [-9 , 9 }; 



with 



\C(re l9 )-C(r)-i9a(r)+9 2 b(r)/2\<\9 3 \/6 sup |/ (3) («)l 

\a\<6 



f (3) (a) = \~ire la C'(re ia )-3ir 2 e 2la C"(re la )-ir 3 e 3la C'"(re' a ) 



< rC'{r) + 3r 2 C"{r) + r s C"'{r) 



(l_ re )7/2 

as r — > 1/e, for some constant re. Hence, if we take 9$ = 9a(r) = o((l — re) 7 / 6 ), then 

A{re lB ) 



sup 

|0|<0o(r) 



-i8a(r)+9^b(r)/2 



A{r) 







(65) 



as r — > 1/e. Thus H2 holds for such values of #o- We now take 

o (r) = (l-re) 6 / 5 , 

and want to prove that H3 also holds. 

Recall that C(z) is analytic on C \ [1/e, 00), and let us isolate in C(z) the part that diverges 
at z = 1/e: 

CW-^ + ilog^+OCl) (66) 
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where c = e - 5/4/V2 > 0. It follows that 

B{z) := C{z) - - \ log 

Vl — ze 4 1 — ze 

is uniformly bounded on {\z\ < 1/e}. Hence, writing z = re l ° , we have 



sup 



A(z) 



< M sup 

|0|e[0o,O 



1 — re 


1/4 


1 — ze 





exp 



Vl — ze vl — re 



A(r) 

for some constant M. 

For any z of modulus r < 1/e, we have |1 — re| < |1 — ze\, and we can bound the first factor 
above by 1. Also, it is not hard to prove that 3? (l/Vl — ze) is a decreasing function of 9 £ (0, it). 
Hence, denoting z — re l9 ° : 
A(z) 



sup 

\e\e[6 ,< 



A(r) 



'b{i 



< M exp !ft 



But as r — ¥ 1/e, the choice (65 1 of 9q implies that 



y/l - zoej Vl - re 
o(l). 



<b{r). 



3c 



10 



VI - z e J y/T^Fe~ 8(1 -re) 1 / 
Condition H 3 now follows, using the estimate (64) of b(r). 

Let us now consider the series Ai{z) := C ( - l+1 \z)A(z), for i > 1. It is easy to prove by 
induction on i that for i > 1, 



C w (z) = 



(2i)!(e-5/4) £ 



1 



(67) 



4*V2i!(l-ze) l + 1 /2 ' ~ \(1 - ze)* 

This can be proved either from the expression of C(z) given in Proposition [30j or by starting 
from the singular expansion (66) of C(z) and applying [191 Thm. VI. 8, p. 419]. 

Recall the behaviour (64) of the functions a(r) and b(r) associated with A(z). It follows, with 
obvious notation, that asr-> 1/e, 

1 



ai {r) = o(r) + i' " = o(r) + O 



c(*+!)(r) 



1 — re 



and 



. , C( l+2 )(r) 2 C( l+3 )(r) 2 



C(''+ 2 )(r) 



= 6(r) + O 



1 

(1 - re)' 



(68) 
(69) 



both tend to infinity. Thus H 1 holds. 

Let us now prove that Ai (z) satisfies H2 with the same value of 0q as for A(z) (that is, 
6» = (1 - re) 6 / 5 ). Thanks to (|68}|69j) we have, for |0| < 6 and uniformly in 9, 

e -i0a z (r)+e 2 b z (r)/2 = e ~i9a(r)+e 2 b(r) /2 ^ + _ re )l/5)) _ 



Now using (|67|, we obtain, denoting z = re , 

-i-3/2 

(l + 0((1 - re) 1 / 2 )) = 1 + 0((1 - re) 1/5 ). 



c (»+i)(. 



Hence 



C( l +!)(r) 
A ( (z) 



1 — ze 
1 — re 



-iea i (r)+6» 2 /2fc i (r) 



^ e -ifla(r)+flV36(r) (l + ((1 - re) 1 / 5 )) 



Mr) A(r) 
and Condition H2 holds for Ai since it holds for A. 

Finally, since C(z) has non-negative coefficients, we have |C^ +1 '(z)| < C^ l+1 \r) for z = re 10 . 
Thus the fact that Ai(z) satisfies H3 follows from the fact that A(z) satisfies H3, together with 
bi(r) ~ b(r). . 
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Appendix B. Forests with a stable set: a parametrized converging example 

We study in this appendix a family A of graphs that slightly differs, in nature, from the 
examples studied so far: our graphs are now coloured forests, with black and white vertices, 
with the condition that two black vertices are never adjacent. That is, black vertices form a 
stable set. Moreover, we do not only consider the uniform distribution on the set A n of graphs of 
size n, but, more generally, we take each graph G of A n with probability proportional to y b ^ G \ 
where y > and b(G) is the number of black vertices in G. This leads us to introduce bivariate 
generating functions, like 

A (*>v) = E Qy b{G) = E z ^(y), 

and similarly for the class C of trees of A. The main interest of this class lies in the fact that 
it can be made "more or less addable" (in the sense of [30 ) by tuning the parameter y. As we 
shall see, the series C(z,y) converges at its radius p (which now depends on y), and we shall 
prove that most results of Section [3] still hold, except that the probability that a graph of size n 
is connected tends to a limit which, depending on y, can be arbitrarily small. 

Proposition 41 (The number of forests with a stable set). The generating function T Q (z) 
and T m (z) counting rooted trees with a white root and a black root respectively satisfy 

T (z) = ze T »W +T 'W, 

T.{z) = yze T ° {z) . 

The number of (unrooted) trees of size n with j black vertices is 



and the associated generating function is 

C ^y^ = E ~",\ u ' r "- = T °W + W - To(z)T.{z) - l -T {zf. (70) 
The generating function of forests of A is A(z,y) — e c ^ z,v K 

Proof. The system of equations defining T (z) and T,(z) is classically obtained by decomposing 
the trees at their roots. The bivariate Lagrange inversion formula (see, e.g., [2D]), applied to this 
system, gives the numbers t° n a and t* j of rooted trees of size n having j black vertices, and a 
white root or a black root, respectively: 

Hence the number of rooted trees of size n having j black vertices is 

t n>j : M>, .,)■ '. 

Divising by n (the number of possible rootings) gives the above expression of c nj . This expres- 
sion can also be obtained bijectively by adapting several known proofs dealing with properly 
bicoloured trees. Finally, it is also a special case of [101 Thm. 3]. 

Let us now prove the expression of C(z, y). Each unrooted tree of size n is counted n times 
in T(z, y) := T {z) + T,(z). Such a tree has n — 1 edges, so we need to substract from T(z, y) 
the generating function of trees with a marked edge to obtain C{z, y). But a tree with a marked 
edge decomposes into two rooted trees, not both rooted at a black vertex. This gives the term 



-T T. - \Tl in the expression (70 1 of C(z,y). 
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Proposition 42 (Asymptotic enumeration). For y > 0, let j3 be the only real of [1, +oo) 
such that 

y = P{fi-l)eP- 1 . (71) 
Then C(z,y) and A(z,y) have radius of convergence 

P = 

with j3 := 1//3. Moreover, 

Zy/TT 



with 7 = ft\J p2+^ i ■ The probability that Q n is connected is 

c «^) MA, ^ t y ; +h ri \ 2/3 3 -4/3 2 + 4/3-1 

—^-y ~ 1/A(p, y) = exp (-a) with a = C{p,y) = — 2 . (72) 

As y increases from to oo, the parameter (3 increases from 1 to oo, and the above limiting 
probability decreases from l/y/e to 0. 

Proof. Let us still denote T = T(z, y) = T (z) + T,(z). This series can be described by a single 
equation: 

T = G{z,T) 

where 

G(z, w) = ze w + yz exp (ze w ) 
is an entire function of z and w. It is not hard to check that T satisfies the smooth implicit 
function schema of |19l Def. VII. 4, p. 467]. The characteristic system G(p,s) = s, G' w (p,s) = 1 
holds for p — /3e 1 ~' 9 ~^ and s — (3 + (3 — 1 once y is parametrized by (71 1. Theorem VII. 3 of 
p. 468] thus applies: the series T(z, y) has a square root singularity at p, and one can compute 
its singular expansion at this point to an arbitrary ordeij^J 

T(Z, y)=p + P-\- 1 y/1 - z p - ^ ^ ^ 

3(/3 2 +/3-l) 2 

7 (_9 + 77/ 3- 97/3 2 -10/3 3 -7/3 4 + 25/3 5 +30/3 6 + 5/3 7 )(l-z/p) 3/2 

3 r- O ((1 - Zip) ) 

36 (^ 2 +/3- i f V ; 

where 7 is given in the proposition. We then use T a = ze T and T, = T — T Q to obtain expansions 
of these two series at p. Finally, getting back to (70 1, we obtain for C(z,y) an expansion valid 
in a A-domain: 

C(z, y) = a-(P + f3- 1)(1 - z/p) + ^(1 - z/pf 2 + 0((1 - z/p) 2 ) (73) 



= p C(z,y). 



where a is given by (72 1. Our final step is the expansion of A(z, y) = e 

A(z, y) =e"(l-(P + l3- 1)(1 - z/p) + y (1 ~ ^/p) 3/2 + 0((1 - z/p) 2 )) . (74) 

The estimates of c n (y) and a n (y) now follow from singularity analysis. a 

Proposition 43 (Number of components — forests with a stable set). Fix y > 0, and 

define (3 and a as in Proposition \42\ The random variable N n — 1 converges in law to a Poisson 
distribution of parameter C(p,y) = a (see 

Proof. We start from 

[z n }c(z, y y 



P(N n = i) 



i\[z n }A(z,y) 



and apply singularity analysis using the expansions (T73J) and (74 1 of C(z,y) and A(z,y) 



3 This is only done at first order in 1191 Thm. VII. 3], but the possibility of pushing the expansion further is 
mentioned at the top of p. 469. 
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Proposition 44 (Size of the root component — forests with a stable set). The random 
variable n — S n converges to a discrete limit distribution X given by 



where p, and A(z, y) are given in Proposition 41 . The same holds for the size L n of the 
largest component. 

Proof. The two ingredients used in the proof of Proposition [5] namely smoothness of a n {y)/n\ 
and convergence of c n (y)/a n (y), still hold here. ■ 
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